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the charm quark propagating as an overlap fermion at fixed lattice spacing. The dispersion relation 
and mass splittings are evaluated. 
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D s spectrum in 2+1 DWF with overlap charm 



1. Introduction 

The discoveries of new resonances D s j by the B factory experiments ^ and CLEO [||] have 
provoked much interest in heavy-light systems in general and in the D s mesons in particular. The 
mass splittings can be understood in terms of heavy quark and chiral symmetry |[] . 

In the double limit of heavy quark and chiral symmetry, the two heavy-light multiplets, {0~ , l - } 
and {0 + , 1 + }, are degenerate. Then chiral symmetry breaking causes splitting between parity part- 
ners, such that the 1 + — 1 and + — 0~ are equal. Experimentally, the splittings, shown in Table 
|l[ are remarkably close. The hyperfine splitting can also be understood in terms of heavy quark 
symmetry breaking effects. 

Many previous lattice calculations ^ [7|, pi |9|, [10|, [IT]] tried to reproduce the features of these 
heavy mesons, most of them considering a static or non-relativistic heavy charm quark, with the 



exception of [g|] which uses the Fermilab approach and [ |11| ] which describes the charm quark as a 



domain wall fermion. All these works are in the quenched approximation. 



In this work the charm quark is described by an overlap Q12Q formalism, while the light strange 
quark is a domain wall fermion, DWF [|l3|]. 

0+-0 1+-1- 1-0 



349.1(4) 346.9(1.0) 143.8(4) 
Table 1: Experimental values in MeV for the different mass splittings from 



2. Numerical details 

The gauge ensembles used for our calculations are the 2+1 flavour dynamical DWF ensem- 



bles from RBC-UKQCD collaboration [14]. They were generated with the renormalized group 
improved Iwasaki gauge action at j8 = 2.13. The lattice volume is 16 3 x 32, with the fifth di- 
mension L s = 16 and the domain wall height aM^ = 1.8. Three sea quark masses are considered, 
am sea = 0.01,0.02,0.03, and the strange quark mass is fixed at am s = 0.04 [|15|]. The correlators 
were measured with sources on multiple time planes, in order to improve our statistics. Details of 
the three ensembles used are listed in Table ^| 

For the overlap charm [[[(]] quark mass, two values are chosen, am c ~ 0.72, 0.9. Correspond- 
ingly, we have two heavy-light mesons, indicated as HI, the lighter, and H2, the heavier, for each 
sea quark mass. Recall the expression of the massive overlap operator: 

aD 0V = p(l + n) + p(l - n)y 5 sgn(y 5 (aD w - p)), (2.1) 

where /I = and p is any mass parameter that can be added to Dw without affecting the con- 
tinuum limit: here it was chosen equal to 1.3 looking at the heavy-heavy pseudoscalar. The over- 
lap operator was used to invert on hyp-smeared DWF gauge configurations for mass parameter 
jU ~ 0.277, 0.346, corresponding to the two charm mass values above. 
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am sea N tm j(sep) No. of Origins p 2 m N cfgs 



0.01 


500-4000(50) 


4 


4 


282 


0.02 


1000-4025(50) 


2 


4 


122 


0.03 


1000-4000(50) 


2 


2 


121 



Table 2: Datasets used. 



3. Analysis 

In Figure [j] we show typical effective masses for the low-lying J p states of the four channels 
we are interested in. The left plot is for the heavy-light meson containing the lighter charm quark, 
the right plot for the heavier one. For the pseudoscalar and vector channels, similarly reasonable 
plateau are found for higher momenta. Once computed the meson masses at different lattice mo- 
menta, 

p L a = ^, nGZ + (3.1) 
La 

we fit them to the dispersion relation. 

The dispersion relation is defined such that the iff{m 2 a 2 ) error is reflected in the deviation of c, 
the effective speed of light, from unity. We fit the energies to a quadratic expression as in eq. (3.2), 
(3.4) , as explained below. 



The plot in Fig. ^| shows the dispersion relation for the pseudoscalar lightest meson, i.e. HI, 
in the am sea = 0.01 case with 5 momenta. The value of the speed of light obtained from the 
fit to eq. (3A), c = 0.897(12), is higher than one might expect from One of the methods 



trying to overcome the problems with heavy quarks (i.e. aniQ ~ 1) is the Fermilab or Relativistic 



Heavy Quark approach [17]. It gives us an alternative interpretation of the dispersion relation: the 
basic idea is considering the expansion of the energy-momentum relation in powers of (lattice) 
momentum pa, 

{Eaf = (M ia ) 2 + Ml(pa) 2 +K(pa) 4 + ... (3.2) 

M2 

where M\ is the rest mass, M\ =E(0), and Mi is the so-called kinetic mass, = ( j4)p=o- The 
relativistic mass shell will have niQ =M\= M2, and the expression above becomes 

{Eaf = (Mia) 2 + (c 2 = \){paf . (3-3) 

In practice, at our non-relativistic mass, we can not truncate the expansion at p 2 , but we have 
to consider higher order terms, i.e. including 8Ei at , 

{Eaf = {m Q af + c 2 {paf + 8E lat . (3.4) 

It has been observed that the rest mass of non-relativistic particles decouples from the inter- 
esting dynamics. 



'Only three momenta are available for the am sea = 0.03 ensemble, i.e. there are only three points in the plot of 
energies versus momenta, so we found the linear fit with pa = 2sin(^y / n/Lfl) (following [|lrj]) be the best one in the 
0.03 case. 
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The suggestion from the Fermilab approach [Q] is then considering M2 instead of M\ and 
tuning the couplings in the lagrangian so that M2 takes the physical value. In this preliminary 
analysis we consider both Mi and M2 and look at the dependence of the mass splittings on them. 
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Figure 1: Effective mass for the low-lying heavy-light mesons with charm mass am c ~ 0.72 (left) and 
am c ~ 0.9 (right), for the am v(M = 0.01 ensemble. 




Figure 2: Dispersion relation plot for the pseudoscalar channel of the lighest heavy-light meson considered, 
in the m = 0.01 case. 



4. Results 

First of all let's clarify the notation used for the mass splittings considered: AH = 1 — is 
the hyperfine splitting, AS = + — 0~ and AV = 1 + — 1~ are the scalar and vector parity splitting 
respectively. The values of these splittings obtained with all three ensembles for both our heavy- 
light mesons, HI and H2, are summarized in Table |3[ The same splittings values are plotted versus 
am sea in Figure ||. We can notice a very small dependence on the sea quark masses. 

The plots in Figure ^ summarize our main results. The two plots on the left show our splitting 
values versus 1 /Mps, where Mp$ is equal to Mi in the upper left panel and to M2 in the lower left 
one. The two plots on the right show the ratio of vector parity over scalar parity, 1 + — l~/0 + — 0~ 
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m sea Meson aAH aAS aAV 



0.01 


HI 


.109(5) 


.265(20) 


.255(31) 




H2 


.108(8) 


.260(25) 


.240(52) 


0.02 


HI 


.121(8) 


.202(53) 


.173(65) 




H2 


.125(12) 


.274(32) 


.228(37) 


0.03 


HI 


.131(6) 


.310(67) 


.232(46) 




H2 


.110(14) 


.267(43) 


.256(34) 



Table 3: Mass splittings value obtained for each ensemble and for both our heavy-light mesons, HI and H2. 



Figure 3: Mass splitting values listed in Table 3 versus the am a 



♦ AH for HI 

O AH for H2 

♦ AS for HI 
O AS for H2 

♦ AVforHl 
© AV for H2 



0.02 

am 



versus l/Mp$, as before. The horizontal line represents the experimental value. Results from all 
three m sea values are shown. In all plots, the vertical line represents our estimate of the physical 
D s meson, using aT x = 1.60(3) Gev [15]. We can see that the effect of using M2 instead of Mi is a 
shift in the x axis, as we expect from eq. (3.2) and (^4): the difference between the two masses is 
entirely a lattice artefact (eq. 3.4). 

What we can see from these four plots is in any case no heavy quark mass dependence of the 
mass splittings. 



5. Conclusions 

In the first stage of our study of D s meson on 2+1 DWF QCD, with the charm quark as an 
overlap fermion, we found clear signals for all the four channels we were interested in. A very 
little dependence of the splittings on the sea quark mass is observed. For the dispersion relation 
analysis, both M\ and M2 were considered: as expected, no heavy quark mass dependance in the 
splittings is observed. The ratio of two parity splittings obtained is close to the experimental value 
within statistical errors, as shown in Figure ||. The am sea = 0.03 ensemble data don't always follow 
the trend of the other two: investigations are in progress. In order to reduce our large error bars, we 
need more statistics. A quenched calculation is also in progress. 
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AH for H 1 
AH for H2 
AS for HI 
AS for H2 
AV for HI 
AV for H2 



0.82 0.84 
1 / (aMJ 
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AH for H2 
AS for HI 
AS for H2 
AV for HI 
AV for H2 



0.4S 0.5 0.55 0.6 0.65 0.7 0.75 0.8 0.85 0.9 
U(aM ) 
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0.8 0.82 0.84 
1 / (aM t ) 



0.86 0.88 0.9 




0.45 0.5 0.55 0.6 0.65 0.7 0.75 0.8 0.85 0.9 



Figure 4: On the left the splitting values versus 1 /Mp$ are plotted, with Mp$ equal to M\ in the upper left 
panel and to M2 in the lower left one. On the right the ratio 1 + — l~/0 + — 0~ obtained is plotted versus 
1/Mps, with Mps as before. The horizontal line corresponds to the experimental value. In all plots, the 
vertical line represents our estimate of the physical D s meson, using aT 1 = 1.60(3) Gev jl5|]. 
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Figure 5: Plot of the 1 + — 1 /0 + — ratio obtained for the three ensembles versus the hyperfine splitting 
over Mi . The experimental value is also plotted. 
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